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Abstract 



The common interpretation of HBT data measured at top SPS energies leads to apparent source 
lifetimes of 6-8 fm/c and emission duration of approximately 2-3 fm/c. We investigate a scenario 
with continuous pion emission from a long-lived (~ 16 fm/c) thermalized source in order to show 
that it is not excluded by the data. Starting from a description of the source's spacetime expansion 
based on gross thermodynamical properties of hot matter (which is able to describe a number 
of experimental observables) , we introduce the pion emission function with a contribution from 
continuous emission during the source's lifetime and another contribution from the final breakup 
and proceed by calculating the HBT parameters R ou t and Rside- The results are compared with 
experimental data measured at SPS for 158 AGeV central Pb-Pb collisions. We achieve agreement 
with data, provided that some minor modifications of the fireball evolution scenario are made and find 
that the parameter R ou t is not sensitive to the fireball lifetime, but only to the duration of the final 
breakup, in spite of the fact that emission takes place throughout the whole lifetime. We explicitly 
demonstrate that those findings do not alter previous results obtained within this framework. 



1 Introduction 

One of the main goals in the study of ultrarelativistic heavy-ion collisions (URHIC) is to understand the 
conditions realized in the center of such a collision. If the created system is sufficiently thermalized, one 
can in principle link measurements with information coming from lattice simulations about the equation 
of state (EoS) of hot QCD matter (see e.g. ^]). However, in order to make use of the EoS, some 
information on the spatial expansion pattern of the matter is extremely helpful. 

HBT correlation measurements can provide at least part of this information (see [21 121 f° r a review). 
However, such measurements can only reveal the correlation lengths of the source, which result from a 
complex interplay of temperature, flow pattern and true geometrical shape. Commonly, simple (often 
Gaussian) parameterizations of the emission source are used to interpret the experimental information 
and extract the geometrical (Gaussian) radius of the source Rq, its transverse expansion velocity v±, 
from the longitudinal correlation radius Ri ong the source lifetime r/ and from R ou tl Rside its emission 
duration At. In 0] (discussing HBT measurements of Pb-Pb collisions at SPS), these parameters are 
found to be Rq ss 6 fm, ry « 6 — 8 fm/c, v± w 0.5c and At rj 2 — 3 fm/c for an assumed temperature 
at breakup of Tf = 120 MeV. These numbers point to a rather rapidly expanding and decaying system. 

However, there is a subtle problem with this interpretation: If the system expands from a root mean 
sqare radius i? rms of about 4.5 fm (as found by a calculation of the overlap of the colliding nuclei) to the 
measured R rm s of about 8.5 fm (assuming Rq = 6 fm), while at the same time the fireball front (for a 
Gaussian distribution, take e.g. R rm s for the front) expands to a velocity of ~ 0.5c, we may expect the 
following equations to approximate the problem |1U|: 



Rrmsirf^) Rrmsify 




and u_l(t/0 = a ■ Tf 



(1) 
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Solving for a and Tf we find Tf « 15 fm/c, clearly incompatible with the lifetime extracted from Ri ong . 
What is wrong here? 

In [TJ]|, it has been pointed out that the relation [5] 

fliong = Tf(T f /mt) 1/2 (2) 

(with transverse mass m t = ypf-h m 2 and pt the transverse momentum of a particle with mass m) 
linking the measured Ri on g and t/ is based on a backwards extrapolation of the measured freeze-out 
state, assuming that the longitudinal expansion takes place with the same velocity at all times before 
Tf. However, this relation ceases to hold if the system undergoes accelerated longitudinal expansion also 
(in this case, the lifetime can be significantly longer than extracted from @), so the estimate done from 
the transverse dynamics is probably more reliable. 

This, however, immediately gives rise to the following question: If a thermal source expands for 0(15) 
fm/c, how can it be that its apparent emission duration is only 0(2) fm/c? In order to answer this ques- 
tion, one needs to specify the space-time evolution of the source beyond the schematic parameterization 
used in 0]. We will employ a description of this evolution based on bulk thermodynamic properties of 
the fireball to investigate this question. 

The paper is organized as follows: After an introduction to the fireball expansion model, we discuss in 
greater detail why longitudinal acceleration helps to match the timescales extracted from Ri ong and R S ide ■ 
Then we specify how this model is used to calculate the relevant HBT parameters. We subsequently 
compare the results to the measured data and explicitly demonstrate the differences to the standard 
parameterizations. Then we discuss the modifications to our scenario which are required to obtain good 
agreement with the data. Finally we investigate the physics behind the observed long fireball lifetime 
and the small difference between R out and R S ide and discuss the impact of the modifications introduced 
before to previous results obtained within the same framework. 

2 Fireball evolution and hadron emission 

In this section, we briefly outline the general framework of the fireball evolution model which we use in 
the following as a starting point for the calculation of hadron emission and ultimately the calculation of 
HBT radius parameters. 

2.1 Expansion and flow 

Our fundamental assumption is to treat the fireball matter as thermalized from an initial proper time 
scale To until breakup time Tf. For simplicity, we assume a spatially homogeneous distribution of 
matter. Since some volume elements move with relativistic velocities, it is sensible to choose volumes 
corresponding to a given proper time r for the calculation of thermodynamics, hence the thermodynamic 
parameters temperature T, entropy density s, pressure p, chemical potentials ^ and energy density e 
become functions of r only for such a system. In the following, we refer to r as the time measured in a 
frame co-moving with a given volume element. 

In order to make use of the information coming from lattice QCD calculations, we proceed by calculating 
the thermodynamical response to a volume expansion that is parametrized in such a way as to reproduce 
the experimental information about the flow pattern and HBT correlations as closely as possible. 

As a further simplification, we assume the volume to be cylindrically symmetric around the beam (z)- 
axis. Thus, the volume is characterized by the longitudinal extension L(t) and the transverse radius 
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i?(r) and we find 



V(t) = ttL{t)R 2 (t). 



(3) 



In order to account for collective flow effects, we boost individual volume elements according to a 
position-dependent velocity field. For the transverse flow, we make the ansatz 

VT^Tj^r/Rr^iT^iT) (4) 

where -R rms (r) denotes the root mean square radius of the fireball at r and f?™ s (i") the transverse 
rapidity at i? rms . 



For the longitudinal dynamics, we start with the experimentally measured width of the rapidity interval 

7? 



of observed hadrons 2ny ont at breakup. From this, we compute the longitudinal velocity of the fireball 



front at kinetic freeze-out vj ront . We do not require the initial expansion velocity VQ ront to coincide 

with vj ront but instead allow for a longitudinally accelerated expansion. This implies that during the 
evolution r\ = r\ s is not valid (with n — atanh v and rj s the spacetime rapidity i} a = 1/2 ln(i + z)/(t — z)) 
in contrast to the non-accelerated case. 

The requirement that the acceleration should be a function of r and that the system stays spatially 
homogeneous for all t determines the velocity field uniquely once the motion of the front is specified. 
We solve the resulting equations numerically [S] and find that for not too large rapidities r\ < 4 and 
accelerations volume elements approximately fall on curves const. — \/t 2 — z 2 and that the flow pattern 
can be approximated by a linear relationship between rapidity r\ and spacetime rapidity r\ s as rj{rj s ) = ^r/ s 
where £ = r]^ ront /rjl ront and rj^ ront is the rapidity of the cylinder front. In this case, the longitudinal 
extension can be found calculating the invariant volume V — J da^u^ as 

sinh ((C - l)n front (T)) 
L(t) « 2t U V j / * ^ (5) 

(C-i) 

with r]( ront (T) the spacetime rapidity of the cylinder front. This is an approximate generalization of 
the boost-invariant relation L(t) — 2r\$ Tont T which can be derived for non-accelerated motion. A more 
detailed description of the expansion can be found in 6 . 



2.2 Parameters of the expansion 

In order to proceed, we have to specify the longitudinal acceleration a z (r) at the fireball front (which 
in turn is used to calculate r)( ront (T) numerically), the initial front velocity VQ ront and the expansion of 
the radius R(t) in proper time. 

We make the ansatz 

P( T ) ( c\ 

a z = c z - —— (6) 
e(rj 

which allows a soft point in the EoS where the ratio p/e gets small to influence the acceleration pattern. 
c z and VQ Tont are model parameters governing the longitudinal expansion and fit to data. 

Since typically longitudinal expansion is characterized by larger velocities than transverse expansion, i.e. 
v front s^, w / ront j we treat the radial expansion non-relativistically. We assume that the radius of the 
cylinder can be written as 

R(r) =R + c T [ T dr> f dr"^Q- (7) 

Jt Jto £ ( t ) 

The initial radius Rq is taken from overlap calculations. This leaves a parameter ct determining the 
strength of transverse acceleration which is also fit to data. The final parameter characterizing the 
expansion is its endpoint given by t/, the breakup proper time of the system. 
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2.3 Thermodynamics 



We assume that entropy is conserved throughout the thermalized expansion phase. Therefore, we start 
by fixing the entropy per baryon from the number of produced particles per unit rapidity (see e.g. 0). 
Calculating the number of participant baryons (see [8]) we find the total entropy Sq. The entropy density 
at a given proper time is then determined by s = Sq/V(t). 

We describe the EoS in the partonic phase by a quasiparticle interpretation of lattice data which has 
been shown to reproduce lattice results both at vanishing baryochemical potential [Ib and finite fis |H| 
(see these references for details of the model). 

In the hadron gas phase, we calculate thermodynamic properties at kinetic decoupling where interactions 
cease to be important. Here, we have reason to expect that an ideal gas will be a good description. 
Determining the EoS at this point, we choose a smooth interpolation from here to the EoS obtained in 
the quasiparticle description. This is described in greater detail in [HI E] • 

With the help of the EoS and s(t), we are now in a position to compute the parameters p(t), e(r), T(r) 
as well. Since the ratio p(r)/e(r) appear in the expansion parametrization, we have to solve the model 
self-consistently. 



2.4 Solving the model 

In order to adjust the model parameters, we compare with data on transverse momentum spectra 
and HBT correlation measurements. In a similar model is fit to a large set of experimental data, 
providing sets of parameters Tf,v±f,Rf, r)j ront at freeze-out. Although the present model uses a different 
(box vs. Gaussian) longitudinal distribution of matter, we use the parameters from this analysis as a 
guideline for the transverse dynamics where this difference should not show up and determine rf™ nt 
separately. Specifically, we use the set bl from ^U] for the transverse dynamics. 

By requiring R(rf) — Rf and v^ ront — v±f we can determine the model parameters ct and Tf. c z 
is fixed by the requirement rj^ ront {Tf) — r\^ ront . The remaining parameter VQ Tont now determines the 
volume (and hence temperature) at freeze-out and can be adjusted such that T(rf) = Tf. 

The model for 5% central 158 AGeV Pb-Pb collisions at SPS is characterized by the following scales: 
Initial long, expansion velocity v^ ront = 0.5c, thermalization time tq = 1 fm/c, initial temperature 
To = 305 MeV, duration of the QGP phase tqgp = 7 fm/c, duration of the hadronic phase Th a d = 9 
fm/c, total lifetime Tf — tq = 16 fm/c, r.m.s radius at freeze-out = 8.55 fm, transverse expansion 

velocity v±f — 0.537c. 



2.5 Initial compression and re-expansion 

We have chosen the above model framework in such a way that it allows to account for the apparent 
mismatch in timcscales derived from Eqs. Q and J2J- This is done by allowing for initial compression 
of the matter in longitudinal direction and subsequent re-expansion. The longitudinal flow pattern of 
matter is schematically indicated in Fig. ^ 

It is apparent from the figure that a measurement of the conditions around the freeze-out hypersurface 
t = Tf (where the bulk of hadronic emission takes place) would not be able to distinguish between an 
accelerated or non-accelerated longitudinal evolution history. This has also been pointed out in |1U|. 
Hence, incorporating accelerated longitudinal expansion we have the freedom to select r/ accoring to 
Eq. Q without being in contradiction with (J2J). (Since we in practice adjust the model to Ri on g and the 
observed multiplicity distribution in rapidity, a detailed discussion of Ri ong in this paper is not fruitful) . 
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Figure 1: A scetch of the accelerated longitudinal expansion pattern of the fireball matter. 

In our model, the parameter controlling the initial amount of compression is VQ Tont = 0.5c. This has 
profound consequences for the initial state: If we set VQ Tont = Vj° nt (and hence use the scenario 
underlying Eq. @), the initial temperature of the fireball drops by ~ 30%. 

This in turn has consequences for the emission of real photons. In |12) it was shown that a scenario very 
similar to the one discussed here can account for the measured amount of photon emission, especially 
in the momentum region between 2 and 2.5 GeV where no significant prompt photon contribution 
is expected. Assuming a scenario without longitudinal acceleration, the yield of direct photons at 2 
GeV drops by a factor 8 for the standard choice of the equilibration time of tq s» 1 fm/c in striking 
disagreement with tha data. This cannot fully be compensated even for a choice of the equilibration 
time as low as 0.1 fm/c. 

In conclusion, the longitudinally accelerated pattern advocated earlier to explain the difference between 
the timescales extracted from and is strongly favoured by the direct photon data. Therefore, we 
will assume in the following that the relevant timescale of the fireball evolution can be inferred from 
transverse dynamics only and focus on the remaining question of the apparent mismatch of the source 
evolution time 0(15) fm/c and its emission duration 0(2) fm/c. In order to investigate this, we have to 
specify the emission function of the fireball. 



2.6 Hadron emission from the hypersurface with timelike normal 

One contribution to the total emission of particles from the fireball comes from the final breakup of the 
thcrmalized system. The corresponding freeze-out hypersurface is characterized by a timelike normal. 
We calculate particle emission from this phase with the help of the emission function 



. | Mt cosh(F — 7] s ) f K ■ u(x) + fj,i 
tcLt / (t — Tf) 2 " 



S\{x, K)d 4 x ^ ^ exp ^ ± 1 ) • G(r)H( Vs )d Vs rdrd<P 



-l 



: exp 



v /2tt(At) 2 V 2 ( At ) 



(8) 



In this expression, the momentum K is parametrized in terms of longitudinal rapidity Y , transverse 
mass Mt = \J Kf + mf and transverse momentum K t The factor exp ^ fc '"^ +M ' ± 1^ corresponds to 
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the thermal distribution of particle species i (with corresponding chemical potential fii and mass rrij) 
where the product k ■ u(x) takes care of evaluating the distribution in the rest frame of matter moving 
with four- velocity u(x) as seen from a particle with momentum K and the +(— ) sign holds for fermions 
(bosons). G(r) and H(rj s ) describe the distribution of matter in radial and longitudinal direction using 
the radius r and spacetime rapidity r) s . Finally, the factor exp (— (r — t/) 2 /(2(At) 2 )) can be seen as 
a smearing of a sudden decoupling occuring at proper time r/ into an extended breakup time period 
characterized by At. 

The single particle spectrum for particle species i is calculated from this emission function as 

K^=J d±xSl(x,p) (9) 
with Ep = \/ p 2 + m 2 . HBT correlation radii are obtained using the common Cartesian parametrization 

C(q, K)-l = exp [-q 2 R 2 0Ut (K) - q 2 s R 2 slde (K) - q?Rl ng (K) - 2q oqi R 2 ol (K)] (10) 

(see e.g. El for an overview and further references) for the correlator. Here, K — \{p\ + P2) is 
the averaged momentum of the correlated pair with individual momenta pi,P2 and q — (pi — P2) the 
momentum difference. The transverse correlation radii R ou t,side follow from the emission function as 

Ride = (f) (11) 



R 2 out = ((2 - [3 J) 2 ) (12) 
with f3± the transverse velocity of the emitted pair, = — (x^) and 

IftrWUa I d i xf(x)S{x,K) 
an average with the emission function. 

In JO], expressions JJjJ, {TJ| and l|12|) have been used to fit the parameters of the shapes G(r) and H{rj), 
the radial and longitudinal expansion velocities at freeze-out 77, the freeze-out temperature Tf and the 
temporal smearing At to experiment, thus providing the basic parameters for the fireball evolution 
described above. As a starting point for the investigations in this paper, we assume a box profile 
G(r) = 9(Rb — r) for the radial distribution of thermalized matter. We use H(rj s ) = 6(ril ront ( T f) — 
Vs)8(v{ ront ( T f) + Vs) throughout this paper. 



2.7 Hadron emission from the hypersurface with spacelike normal 

In addition to particle emission at the final breakup of the fireball described by JSJ) , there is also contin- 
uous emission of particles throughout the whole lifetime of the system. Such emission is characterized 
by a freeze-out hypersurface with spacelike normal. 

Since the correlation radius R ou t probes not only the geometrical radius of the source, but also the 
duration of particle emission (see Eq. I|12|ll. it may well be that such a contribution significantly change 
the ratio Rout I 'Rside (which is known to be small experimentally) as compared to the use of only the 
emission function (JHJ. 
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In order to investigate this question, we consider in addition to © also the emission function 



Sl(x, K)d 4 x = (K± cos(0) - v±(t)M t cosh(F - r) B )) ■ 6 (K x coaU) - v ± {t)M t cosh(F - j] s )) 

exp ( K-^+V? ± ^ drhrdrd(j)TdT S(r ~ R B (r))d(r, s + »^ roBt (r)) 



x exp 



x e(r/ s ront (T) - V s)e(T - n )0{Tf - T) 



(14) 



which describes the emission of particles with transverse velocity K±/Ex cos(0) from the fireball surface 
moving at v± (with (j) the angle between surface and the momentum vector of the emitted particle) at 
r = Rb(t) between the space-time rapidity extent of the fireball from +77/ (t) to — ?7/(t) between initial 

time Ti and freeze-out time r/. The factor 8 {^g^ cos(0) — vj_(t)^ ensures that emission only takes place 
for particles moving faster than the emission surface and no negative contribution of 'backward emission' 
are counted. 

We assume that the expanding fireball can be described as a hot 'core' and a (thin) boundary region from 
which the emission takes place. This region is assumed to be characterized by an average temperature 
Tb- We do not require Tb = Tf, but leave this as a free parameter, however we do require that the 
boundary is characterized by the same pion density as the whole fireball at breakup, which fixes nf for 
given Tb- 

To a first approximation, we describe the total particle emission by the sum of both emission throughout 
lifetime and final breakup as S l (x, k) = S\{x, k) + S l s (x, k). 



3 Calculation of HBT correlation parameters 

3.1 Rside an d the freeze-out geometry 

The experimental results for the correlation radius R S ide are often fit by a function of the form 

Rside = — , ° (15) 
^Jl + m t rf ±f /T f 

This expression comes from an approach where only the final breakup of the fireball is taken into account, 
the matter distribution G(r) is assumed to be Gaussian with width Rq, and Boltzmann statistics instead 
of Bosc-Einstcin statistics is used. 

In contrast, in this paper we study a model with a long duration of the emission (throghout the complete 
fireball lifetime of ~ 17 fm/c), and a final breakup, which is characterized rather by a box-shaped trans- 
verse distribution . In addition, we use the Bose-Einstein (Fermi-Dirac) distribution where appropriate. 

Evidently, the analytic expression (|15|l is only able to determine the ratio T)\f/T. Assuming Tf = 120 
MeV, the fit yields Rq « 7 fm corresponding to R rm s ~ 9.8 fm and v± = 0.55c for the most central 
collisions at 158 AGeV 4,. (Compare with the parameters in our fireball evolution scenario Tf = 100 
MeV, v ± = 0.537c and R rms = 8.55 fm.) 

For the sake of clarity, we start by presenting the result of evaluating (|ll|l using Sf(x,K) only, i.e. we 
neglect all emission throughout the fireball lifetime. The result is shown in Fig. [21 
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Figure 2: Left panel: CERES data on R S ide for central Pb-Pb collisions at 158 AGeV with statistical 
(small bars) and systematical (large bars) errors (unmerged), in addition to a fit to data based on a 
Gaussian source and the result based on the box profile used in our fireball evolution scenario, considering 
breakup at Tf only. Right panel: Data as before, shown is emission only at fireball breakup (dashed), 
and the results with the full emission function for two different values of the average temperature of the 
emission layer (dotted and solid). 



Clearly, the overall normalization (set by the R rm s radius) appears too small in the model and the falloff 
with transverse pair momentum K t too steep. 

The latter effect can be explained by the observation that the Gaussian distribution allows for a (small) 
component of the total amount of matter to be at large radii and hence to be boosted with large radial 
velocities if a flow profile as in Eq. is chosen, whereas the box profile used in the present scenario 
does not. 

This discrepancy is significantly reduced by using the full emission function S v (x, K) instead of S^{x, K). 
This is due to the fact that emission of fast particles is preferred for a freeze-out hypersurface with 
spacelike normal, leading to an enhancement of the correlation radii at larger K t also (Fig. [3 right 
panel). In addition, we have some freedom to adjust Tb (the average temperature of the emitting region 
in SJ(x,K)) within reasonable limits. The data prefer a rather large value Tb ~ 150 MeV. This should 
not come as a surprise, since the temperature of the emitting layer may well be large in the initial stages 
of fireball evolution. 

Guided by the success of the Gaussian parametrization, we drop the unphysically sharp cutoff of the box 
profile G(r) = 6(Rb ox — r) and use instead a Woods-Saxon distribution (retaining the value of Tb = 150 
MeV in the following) 



G(r) = 



1 



exp 



r-R bo .. 
dws 



(16) 



We find that the data prefer rather small values of dws of the order of 0.8 fm (see Fig. [3J left panel). 
This is reassuring since it gives some justification for the use of a box profile in the calculation of dilepton 
and photon emission. 
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Figure 3: Left panel: CERES data on R S ide for central Pb-Pb collisions at 158 AGeV with statistical 
(small bars) and systematical (large bars) errors (unmerged). Left panel: Calaulated R S ide for various 
values of the Woods-Saxon skin thickness parameter dws- Right panel: Calculated R S ide using dws — 
0.8 fm (dashed), shown is the effect of model-independent correction terms <|17[) (solid). 



Finally, we also calculate the model-independent corrections JS] to (|TT1) 

which primarily affect the low momentum region. This is shown in Fig. [3 right panel. These corrections 
are small, however, they tend to improve agreement with data slightly in the low K t region. 

In summary, including the emission throughout the fireball lifetime and all correction terms, we achieve 
a good description of the experimental data. 



3.2 R out and the freeze-out criterion 

It is commonly deduced from and (|12() that the duration of pion emission can be calculated as 

At 2 = -J {R 2 out - Rlide) ( 18 ) 
Pt 

with (] t the velocity of the emitted pair (3 t — K t /Ex- However, there are in fact two different timescales 
related to the duration of the emission in our model: a) the total lifetime of the fireball of order ~ 17 
fm/c and b) the duration of the final breakup, which is a priori unknown, but can be expected to be of 
order of at most a few fm/c for the whole evolution scenario to make sense. We therefore expect the 
actual difference between R ou t and R S ide in our model to be characterized by a mixture of these two 
scales, weighted with their respective contribution to the total pion emission. 

We start the investigation of R ou t by fixing all parameters at the value used to describe R S ide and 
evaluating l|12|) with different values of At, the final breakup duration. In all following calculations, we 
also include the correction term for R ou t |15| 
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Figure 4: CERES data on R ou t for central Pb-Pb collisions at 158 AGeV with statistical (small bars) 
and systematical (large bars) errors (unmerged). Left panel: Calculated R out for different values of the 
fireball breakup time At. Right panel: R ou t for different momentum-dependent freeze-out scenarios (see 
main text) 

The result is shown in Fig. 0] 

Clearly, the overall magnitude of the calculated R ou t is in agreement with data, provided that a rather 
sudden breakup time of the order of 1-3 fm/c is choosen. This indicates that that the difference 
\j0l{R 2 out — R 2 side ) does not reflect the full fireball lifetime r| but rather is dominantly influenced by At 
only (otherwise one would expect the difference to be significantly larger). We will argue later on that 
this is so because the numerical contribtion from continuous emission throughout the fireball lifetime is 
small compared to the contribution from final breakup. 

The shape of the falloff with K t , however, is not well reproduced for any choice of At. (In fact, the shape 
appears similar to the falloff of R S ide and rather insensitive to variations of At.) This suggests that 
the data (albeit the systematic errors are rather large) indicate the need for a more refined freeze-out 
scenario. 

In |16j . such a more detailed scenario based on the (momentum-dependent) mean free path of pions 
inside the hot matter is described. It is beyond the scope of the present paper to show that this freeze- 
out criterion is able to describe the data, however, we will argue that it leads to the correct qualitative 
behaviour of the calculated R ou t- 

The essential observation made in |16j ist that the freeze-out is different for high momentum pions than 
for low momentum pions: High momentum pions are expected to escape earlier. 

If we want to incorporate this into our calculation, we have to make either At or r/ (or both) dependent 
on Kt, the idea being that the final breakup of the fireball occurs earlier for high momentum pions and 
happens faster. 

In the right panel of Fig. we demonstrate that the dependence of At is able to qualitatively improve 
the agreement with data. To obtain an estimate for the qualitative behaviour, we assume At = 2.0 fm/c 
for Kt — 100 MeV and At = 1.0 fm/c for K t — 700 MeV and interpolate linearly inbetween. In view 
of the large systematic uncertainties, however, we do not push this idea further by trying to aim at a 
perfect description of the data. 

As for the dependence of t/ on K t , we observe that earlier freeze-out leads in general to smaller R ou t at 
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larger K t , since one probes a radially less expanded system. A suitable variation of At might then be 
able to lead to a good description of the data. Again, in view of the large uncertainties we refrain from 
a further examination of this scenario. 



4 Pion emission and fireball thermodynamics 

The main result of the previous section is that R ou t does not seem to reflect the order of magnitude 
of the fireball lifetime Tf, even if continuous emission of pions is taken into account. This is a rather 
surprising result, and in this section we will investigate the relevant emission function S^(x, k) further 
in order to answer the following questions: 

• Why is At and not r/ the dominant timescale in i? OM t? 

• Does the continuous loss of pions (and nucleons) have any effect on the thermodynamics of the 
fireball (and hence invalidate the findings of e.g. [5j)? 

In order to answer these questions, we investigate the transverse velocity distributions in the cm. frame 
of pions and nucleons emitted before Tf (which can be easily calculated by integrating Ijl4(l over the 
fireball lifetime). The reason for choosing this particular distribution is the following: We may well 
imagine a process where a particle is emitted in the initial stage of the expansion. Such particle can be 
emitted with relatively small velocity (as seen in the cm. frame), since the fireball initially does not 
expand in radial direction. If the emission velocity of such a particle is, however, smaller than the radial 
velocity of the fireball at r/, there is a chance that this already emitted particle re-enters the thermalized 
region. The condition reads: 

R B (t E )+v^t = R B (tE + t);t E +t<tf (20) 

with emission time and emission velocity ve measured in the cm. frame. This condition is most 
easily implemented using the transverse velocity distribution of emitted particles. 

The resulting transverse velocity distributions are shown in Fig. [S] 




Figure 5: Left panel: Transverse velocity distribution of total it emission and of recaptured ir . Right 
panel: Same for nucleons. Also given is the integral over the distribution. 

Several important observations can be made: For the emission of pions, the recapturing process is not 
very efficient. This is due to the fact that the mass of the pions is of order of the temperature of 



11 



the emitting layer, hence the transverse velocity distribution is peaked towards large velocities where 
recapture cannot take place. On the other hand, the total number of pions emitted before the final 
breakup of the fireball is small as compared to the total number of observed hr (~ 700) . This explains 
partially the observed insensitivity of R ou t to the fireball lifetime. 

Turning to nucleons, we observe that the recapturing process is somewhat more efficient due to the large 
nucleon mass leading to reduced emission velocities. The total number of emitted nucleons, however, is 
also small as compared to the total number of participants (~ 350). 

Note in addition that most of the emission takes place close to 77 anyway, since the area of the emitting 
surface grows strongly with time as the fireball expands and the total number of emitted particles per 
unit time is proportional to this area. This is enhanced by the fact that the fireball in this scenario 
undergoes accelerated expansion in both in longitudinal and transverse direction, thus the emitting 
surface is for the larger part of the evolution significantly smaller than its final value. 

We can incorporate the early loss of pions schematically into the calculation of dilepton emission (see 
[5]) (the observable which is presumably most affected by this effect) by lowering the value of the pion 
chemical potential ^ during the hadronic evolution stage (in we assumed a linear rise with proper 
time of from at the phase transition to a value consistent with the observed number of pions at 
thermal freeze-out. Now, we choose the endpoint such that the number of pions from both continuous 
emission and final breakup agrees with the observed number. We find /^(t/) w 110 MeV). 

The calculation reveals that the dilepton yield in the low invariant mass region between 300 and 700 
MeV is reduced by ~ 10% by the correction (this is hardly visible in a logatithmic plot, therefore we 
do not show a graph). Above 1 GeV invariant mass, the dilepton yield is dominated by Drell-Yan pairs 
and emission from the hot QGP, consequently this domain remains almost unaltered be the correction. 

In summary, we conclude that it is sufficient for the description of gross thermodynamical properties 
of the expanding matter to approximate the fireball as a closed system which does not emit strongly 
interacting partcles before 77. 

5 Conclusions 

We have investigated the effect of particle emission prior to the fireball breakup time 77 in a framework 
based on gross thermodynamic properties of hot QCD matter, which is nevertheless able to account for 
a large number of experimental observables and discussed the connection to the HBT parameters R ou t 
and R S ide- 

Calculating the emission of pions from the fireball boundary and breakup, we computed R ou t and R S ide 
and compared to measurements by the CERES collaboration. We found that R S id e can be well decribed 
within the present framework. 

Regarding R ou t, we identified two timescales which potentially enter this quantity, the fireball breakup 
time At and the total lifetime 77. However, in the present framework, only At was shown to significantly 
influence R ou t ■ In addition, we found that a momentum- independent freeze-out criterion does not provide 
the best description of the data. We argued qualitatively that a more refined criterion (e.g. as suggested 
in |16p improves the agreement with data. 

Finally, we demonstrated that the scenario discussed here does not invalidate our previous results where 
we did not include the effects of continuous emission of particles from the thermalized matter. 
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